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The intercalation phase diagram of Mg in V2O5 from first principles 
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We have investigated Mg intercalation into orthorhombic V2O5, one of only three cathodes known 
to reversibly intercalate Mg ions. By calculating the gronnd state Mg2;V205 configurations and 
by developing a cluster expansion for the configurational disorder in (5-V2O5, a full temperature- 
composition phase diagram is derived. Our calculations indicate an equilibrium phase separating 
behavior between fully demagnesiated a-V205 and fully magnesiated ^-V205, but also motivate 
the existence of potentially metastable solid solution transformation paths in both phases. We find 
significantly better mobility for Mg in the S polymorph suggesting that better performance can be 
achieved by cycling Mg in the 5 phase. 


I. INTRODUCTION 

A multi-valent (MV) battery chemistry, which pairs a 
non-dendrite forming Mg metal anode with a high volt¬ 
age 3 V) intercalation cathode offers a potentially safe 
and inexpensive high energy density storage system with 
the potential to outperform current Li-ion technology.® A 
change in chemistry leads to new challenges, however, one 
being the design of a cathode that can reversibly interca¬ 
late Mg at a high enough voltage. Orthorhombic V 2 O 5 
is one such material that offers exciting prospects of be¬ 
ing a reversible intercalating cathode for Mg batteries.®® 
The theoretical energy density of a cathode based on 
Mg intercalation into V 2 O 5 is ~ 660 Wh/kg,® which ap¬ 
proaches the practical energy densities of current com¬ 
mercial Li-ion chemistries 700 Wh/kg for LiGoOs!^, 
but the major benefit of switching to a MV chemistry is 
the gain in volumetric energy density arising from the us¬ 
age of a metallic anode (~ 3833 mAh/cm^ for MJ^ com¬ 
pared to ~ 800 mAh/cm^ for Li insertion into graphite.®) 

The orthorhombic V 2 O 5 structure has been well char¬ 
acterized due to its interesting spin ladder characteris¬ 
tics and widely known Li intercalation properties, with 
a reversible capacity of ^ 130 mAh/g and voltage of ~ 
3.3 V vs. Li metal.®*!^ Gonsequently, Li intercalation into 
V 2 O 5 has been the subject of several experiment 
and theoreticaPsHUl studies. Li-V 2 05 undergoes several 
first-order phase transformations during intercalation, 
such as the a ^ e and e —i S between xui = 0 and xli = 1, 
the irreversible i5 — >■ 7 transition at xl; > 1 , and another 
irreversible j —i ui transition at xli > 2.1S1 Sev eral a u- 
thors have investigated Mg-insertion into V20^i^*221l25] 
and to date, V 2 O 5 is one of only three cathode materials 
to have shown reversible intercalation of Mg, the other 
two being the chevrel MoaSj^ and layered MoOs.l^ 

While Li-ion has been investigated extensively for the 
past ~ 25 years, there are significantly fewer studies, the¬ 
oretical or otherwise, of Mg intercalation hosts in the lit¬ 
erature. Pereira-Ramos et showed electrochemical 
intercalation of Mg into V 2 O 5 (at 150 °G and 100 /iA/cm^ 
current density), and Gregory et have reported 

chemical insertion of Mg up to Mgo, 66 V 205 . Novak et 


demonstrated reversible electrochemical insertion of 
Mg in V 2 O 5 at room temperature while also demonstrat¬ 
ing superior capacities (~ 170 mAh/g) using an acetoni¬ 
trile (AN) electrolyte containing water as opposed to dry 
AN. Yu et alW^ showed similar improvements in capac¬ 
ity (~ 158.6 mAh/g) using a H2O -I- Polycarbonate (PC) 
system compared to dry PC. Electrochemical insertion 
of Mg into V 2 O 5 nanopowders and thin films using ac¬ 
tivated carbon as the counter electrode was shown by 
Amatucci et and Gershinsky et respectively, 

and insertion into V 2 O 5 single crystals was reported by 
Shklover et 

Thus far, all reported experimental attempts have be¬ 
gun in the charged state and succeeded in reversibly in¬ 
serting only about half a Mg (xMg ~ 0 . 5 ) per formula unit 
of V 2 O 5 , in contrast to M-VgO^where up to XLi ~ 3 has 
been inserted per When the grain size of 

V 2 O 5 is reduced, e.g., nano pow ders a nd thin films, inser¬ 
tion levels can reach XMg 0.6.^^^^ In addition, in cells 
where a Mg metal anode was u sed rapid capacity fade 
was reported upon cycling.^^®^ Unlike Li intercalation 
systems, anode passivation by the electrolytes is a major 
issue for Mg batteries using a Mg metal anode.l^ Out of 
the two experimen ts tha t have not reported significant 
capacity fade so far, 0^P^ lthe work done by Gershinsky et 
al. is particularly useful to benchmark theoretical mod¬ 
els as the Mg insertion was done at extremely low rates 
( 0.5 /rA/cm^), and therefore corresponds most to equi¬ 
librium conditions. 

Previous theoretical studies of the Mg-V 205 system 
have benchmarked structural parameters, average volt¬ 
ages and the electronic properties of layered V 2 O 5 upon 
Mg insertion.l^SHsI] Wang et alW^ showed an increase in 
the Mg binding energy and Li mobility in single-layered 
V 2 O 5 compared to bnlk V 2 O 5 . GarrasccP^ found that 
while incorporating van der Waals dispersion corrections 
in the calculations improved the agreement of the lattice 
parameters with experiments, it led to an overestimation 
of the voltage. Zhou et calculated the band struc¬ 
tures, average voltages, Mg migration barriers, and the 
a ^ 5 phase transformation barrier in Mg-V 2 05 . While 
reporting higher computed average voltage for Mg-V 205 
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compared to the Li-V2 05 sy stem (in apparent disagree¬ 
ment with experimentl;]^g authors explained the 
slow diffusion of Mg in V2O5 by predicting a facile a —?> <5 
transition coupled with an estimated lower Mg mobility 
in S than 

In the present work, we have explored in detail 
the physics of room temperature Mg intercalation in 
orthorhombic V2O5 using first-principles calculations. 
Compared to Li, Mg insertion is accompanied by twice 
the number of electrons, which means that the properties 
of the Mg intercalation system will be largely dictated 
by how the additional electron localizes on the nearby 
V atoms. To study the combined effects not only of 
inserting a different ion but also a different number of 
electrons on the equilibrium phase behavior, we calcu¬ 
late the Mg-V205 intercalation phase diagram using the 
Cluster expansion-Monte Carlo approach. A similar ap¬ 
proach has b een previously used to study Li-intercalation 
system^22133i derived formally through system¬ 

atic coarse graining of the partition function.^^^ Our cal¬ 
culations focus particularly on Mg intercalation into the 
a and S polymorphs of V2O5, evaluating their respective 
ground state hulls, subsequent voltage curves and activa¬ 
tion barriers for Mg diffusion. We have also constructed 
the temperature-composition phase diagram for Mg in 
the S polymorph. 


II. POLYMORPHS OF V2O5 

The V2O5 structure consists of layers of VO5 pyra¬ 
mids, each of which have 4 V—O bonds that form the 
base of the pyramid and one V =0 (Vanadyl) bond that 
forms the apex. Each layer consists of alternate corner 
and edge sharing pyramids, with an offset in the a-axis 
between the edge-sharing pyramids. The different poly¬ 
morphs of V2 O5 observed experimentally are illustrated 
in Figure[iPwith the a (space group Pmmn), 5 (Cmcm) 
and 7 {Pnma) polymorphs all having orthorhombic sym¬ 
metry. The notation, specific to this work, is a being 
the shortest axis of the lattice ( 3.56 A for a; 3.69 A for 
( 5 ), b being the axis perpendicular to the layers indica¬ 
tive of the layer spacing ( 4.37 A; 9.97 A), and c being 
the longest axis ( 11.51 A; 11.02 A). Pure V2O5 crys¬ 
tallizes in the a phase at 298 K and remains stable at 
higher temperatures,^ while the fully magnesiated phase 
(MgV205) has been found to form in the structure of the 
d polymorph.l^ For simpler visualization, a single slice of 
the a and S polymorphs, corresponding to a depth of c /2 
(illustrated by the dashed blue rectangle in Figure [^) is 
shown in Figure and Figure respectively. The a 
and S polymorphs are very similar when viewed along 
the a-axis or the b-c plane (Figure [^). 

The main difference between the S phase and the a 
phase is a translation of alternating V2O5 layers in the a- 
direction by ‘a/2’ which doubles the ‘6’ lattice parameter 
(as well as the unit cell) of the S phase. The Mg sites in 
both a and S are situated near the middle of the VO5 


pyramids (along a) and between the 2 layers (along 6), 
as illustrated by the orange circles in Figure[^ As a result 
of shifting of layers between the a and 6 phases, the anion 
coordination environment of the Mg sites also changes. 
Considering a Mg—O bond length cutoff of 2.5 A, the Mg 
in the a phase is 8-fold coordinated (4 nearest neighbor 
O atoms and 4 next nearest neighbors, 4 - 1 - 4 ) whereas 
the Mg in the 6 phase is 6-fold coordinated ( 4 - 1 - 2 ). In 
this work, the e phase is a specihc ordering of Mg atoms 
on the a-V205 host at half magnesiation, as shown in 
Figure [J. This intercalant ordering is observed in the 
Li-V205 system,ISl and has intercalant ions at alternate 
sites along the a axis , as illustrated by the absence of 
Mg sites in Figure The VO5 pyramids in the a 

and S phases ‘pucker’ upon Li intercalation as observed 
experimentally by Cava et alW^ For the sake of simplicity 
we define puckering here as the angle ‘0’, as shown in 
Figure [^. As the pyramids pucker with intercalation, 
the angle ‘/)’ decreases. 

In the Li-V2 05 system, at xli > 1 , the host structure 
undergoes an irreversible phase transformation to form 
the 7 phase, in which the VO5 pyramids adopt a different 
orientation compared to a and 5 , as seen in Figure 
and[^.I^In the 7 phase, the VO5 pyramids along the c- 
direction alternate between up and down (denoted by ‘-I-’ 
and ‘—’ in Figure]^; whereas, in a and < 5 , the sequence 
goes as ‘up-up-down-down’. The 7 phase has not yet 
been reported in the Mg-V205 system and hence will 
not be further discussed in this paper. 


III. METHODOLOGY 

To compute the ground state hull and the average open 
circuit voltage curves we use Density Functional Theory 
(DFT) as implemented in VASP with the Perdew-Burke- 
Ernzerhof (PBE) exchange-correlation functional.l^^l® 
The Projector Augmented Wave theorj®! together with 
a well converged energy cutoff of 520 eV is used to de¬ 
scribe the wave functions, which are sampled on a F- 
centered 4 x 4 x 4 fc-point mesh. In order to remove the 
spurious self-interaction of the vanadium d-electrons, a 
Hubbard U correction of 3.1 eV is added to the Gen¬ 
eralized Gradie nt Approximation (GGA) Hamiltonian 
(GGA-k[/)PSl as fitted by Jain et al 0 AW Mg-V205 
structures are fully relaxed within 0.25 meV/f.u. 

To obtain the temperature-composition phase dia¬ 
gram, Grand-canonical Monte Garlo (GMC) simulations 
are performed on a cluster expansion (GE) Hamiltonian. 
The GE is a parameterization of the total energy with 
respect to the occupancy of a predefined topology of 
sites , which in this case are the possible Mg insertion 
sites 1241M147I jg practice the GE is written as a truncated 
summation of the Effective Cluster Interactions (ECIs) 
of the pair, triplet, quadruplet and higher order terms as 
given in Equation 
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FIG. 1 : (Color online) (a) a and (b) 5 polymorphs of orthorhombic V2O5 are shown along the c-axis (shown to a depth of 
c /2 for viewing clarity) and along the (c) a-axis, which compared to the (d) 7 polymorph has a different orientation of VO5 
pyramids as denoted by ‘+’ and ‘ —’ signs along the c-axis. Hollow orange circles correspond to the intercalation sites, the green 
dotted lines show the differences in layer stacking and the dashed blue rectangle in (c) indicates a distance of c/2, (e) illustrates 
the e phase corresponding to a specific ordering of Mg atoms in a-V205 at half magnesiation, where alternate intercalant sites 
are occupied in the a axis as indicated by the orange circles. The schematics here correspond to ‘supercells’ of the respective 
polymorph unit cells. 


E{(j) ='^moya(WcFi) ( 1 ) 

OL 

where the energy, E oi a given configuration of Mg 
ions a is obtained as a summation over all symmetrically 
distinct clusters a. Each term in the sum is a product of 
the multiplicity m, the effective cluster interaction (ECI) 
V for a given a, and the occupation variable ai averaged 


over all clusters (3 that are symmetrically equivalent to 
a in the primitive cell of the given lattice. In this work, 
the CE is performed on the Mg sub-lattice and the vari¬ 
ous configurations correspond to the arrangement of Mg 
(cTi = 1 ) and Vacancies (Va; Ui = — 1 ) on the available Mg 
sites. The Pymatgen library is used to generate the vari¬ 
ous Mg-Va arrangements to be calculated with DET.SsHUl 
The CE is built on the DFT formation energy of 97 dis¬ 
tinct Mg-Va configurations using the compressive sens- 
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ing parad igm a nd optimized through the split-Bregman 
algorithm EZEll mean square error (RMSE) and 

the weighted cross-validation (WCV) score are used to 
judge the ciua. lity and the predictive ability of the fit, 
respectively!^ 

The high temperature phase diagram is then obtained 
with GMC calculations on supercells containing at least 
1728 Mg/Va sites (equivalent to a 12 x 6 x 6 supercell of 
the conventional unit cell) and for a minimum of 100,000 
equilibration steps followed by 200,000 sampling steps.l^ 
Monte Carlo scans are done on a range of chemical po¬ 
tentials at different temperatures, and phase transitions 
are detected by discontinuities in Mg concentration and 
energies. In order to remove numerical hysteresis from 
the Monte Carlo simulations, particularly at low temper¬ 
atures, free energy integration is performecP^ with the 
fully magnesiated and fully demagnesiated phases as ref¬ 
erence states. 

Finally, the activation barriers associated with Mg 
diffusion in V2O5 are calculated with DFT using the 
Nudged Elastic Band method (NEB)P^ and forces con¬ 
verged within 100 meV/A. A minimum distance of 9 A 
is introduced between the diffusing species and nine dis¬ 
tinct images are used to capture the diffusion trajectory. 
As previously indicated by Liu et al.^^ the convergence 
of CCA-I- U NEB calculations is problematic, and hence 
standard CCA is used to compute the Mg diffusion bar¬ 
riers. 


IV. RESULTS 

A. Mg-V 205 Ground State Hull 

Figure shows the ground state hull and average volt¬ 
age curves as a function of Mg concentration in V2O5 
as computed by DFT. The solid blue and red lines in 
Figure indicate the ground state hulls of the a and 
5 polymorphs respectively. All formation energies are 
referenced to the fully magnesiated and fully demagne¬ 
siated end points of the d-phase. The overall equilib¬ 
rium behavior of the system is that of phase separation 
between unintercalated a-V205 and fully intercalated 5 - 
MgiV205 as indicated by the solid maroon line. As can 
be observed, the a phase is stable compared to the <5 
phase at low Mg concentrations up to XMg ~ 0.35 where 
the a and 5 hulls intersect, and the 5 phase is stable at 
higher Mg concentrations. In Figure [^, the dash-dotted 
blue line indicates the end members of the a hull (pure 
q;-V 2 05 and a-MgiV205), and the dashed red line the 
lowest energy configurations computed at intermediate 
Mg concentrations for the 5 phase. 

The a-hull represents the energy trajectory for 
metastable Mg insertion into a-V205 (i.e., without trans¬ 
formation of the host to < 5 ), and it displays a convex 
shape with ground state configurations at Mg concen¬ 
trations of 0.25 and 0 . 5 . The most stable configuration 
at XMg = 0.5 in the a hull is the e phase. In contrast. 



FIG. 2 : (Color online) (a) The ground-state hull of Mg in 
V2O5 considering both a and 5 phases. The formation energy 
per formula unit has been plotted with respect to Mg concen¬ 
tration. (b) The average voltage curves at 0 K for the a and 
5 phases with respect to pure Mg metal, obtained from the 
respective hulls are plotted against the Mg concentration. 

there are no metastable Mg orderings in the 5 phase im¬ 
plying that in the d-phase host the Mg ions will want to 
phase separate into MgV205 and V2O5 domains. Some 
Mg configurations when initialized in the a phase relax to 
the 5 phase as indicated by the green diamond points on 
Figure [^. These structures undergo a shear-like trans¬ 
formation from a to d, which involves V2O5 layers slid¬ 
ing along the a-direction. This mechanical instability 
phenomenon has been observed in our calculations both 
at low Mg concentrations (xMg = 0 . 25 ) and at high Mg 
concentrations (xMg = 0 . 75 ), but never at very low Mg 
concentrations (xMg = 0 . 08 ). 

The Mg insertion voltage will depend on which of the 
possible stable or metastable paths the system follows 
and the voltage for several possible scenarios is shown 
in Figure The equilibrium voltage curve is a single 
plateau at 2.52 V vs. Mg metal, consistent with phase 
separating behavior between a-V205 and (5-MgiV205. 
The voltage for the metastable insertion in the a host 
averages ~ 2.27 V vs. Mg metal for 0 < XMg < 1 and ex¬ 
hibits a steep potential drop of ^ 400 mV at XMg = 0 . 5 , 
corresponding to the e ordering. Metastable Mg insertion 
in 6 occurs on a single plateau at 2.56 V vs. Mg metal, 
consistent with phase separation between (5-MgoV205 
and fully intercalated d-MgiV2 05. The average voltage 
of the a phase best agre es with the experimental average 
voltage of - 2.3 V.E^ 


B. Puckering and Layer spacing 

The VO5 pyramids in both a and d-V205 pucker upon 
Mg intercalation, quantified by the angle (j) shown in Fig- 




















5 


ure [^. We find that (j) decreases (corresponding to in¬ 
creased puckering) with increasing Mg concentration, re¬ 
sulting in the formation of ripples in the layers. Current 
calculations show a decrease from (j) ~ 76 ° at XMg = 0 
(which corresponds to flat layers) to (j) ^ 56 ° at XMg = 1 
in the a phase and a decrease from </> ~ 68° at XMg = 0 
to (j) 54 ° at XMg = 1 in 6-Y2O5. 

Figure shows the variation of the V2O5 layer spacing 
(seen in Figure[^ and Figure[^) as a function of Mg con¬ 
centration in both the a (blue) and the <5 (red) phases. 
In other layered materials, van der Waals interactions 
are known to cause layer binding in the deintercalated 
limit , 1 ^ whi ch is not well described by standard DFT 
calculations! ^° l ^^ l Therefore, in order to obtain a better 
estimate of the layer spacing values, additional cal cula- 
tions are performed using the vdW-DF 2 functionalj ^^ * ^^ ! 
which includes the van der Waals interactions in addi¬ 
tion to the Hubbard -I- U Hamiltonian (for removing self¬ 
interaction errors). 



FIG. 3 : (Color online) Variation of layer spacing with 

Mg concentration in both a and <5 phases. The experimen¬ 
tal data points correspond to the pure a-NiOh-, intercalated 
Mgo.2V205 and pure (5-MgiV205. 

The layer spacing values in Figure are taken from 
the relaxed ground states for a and S in Figure [^. 
The blue circles and red squares are obtained from PBE 
(+[/) calculations, while the blue and red triangles are 
calculated with vdW-DF 2 (+[/)■ The experimental val¬ 
ues listed (green diamonds) correspond to pure a-V205,l^ 
Mgo,2V205 rraorted by Pereira-Ramos et al^ and pure 
i5-MgiV205.^As expected, the PBE and vdW-DF 2 layer 
spacing values differ at complete demagnesiation (^ 0.3 
A) but remain similar at all other Mg concentrations, 
where the layer spacing is determined by the electrostat¬ 
ics and short range repulsion. 

With increasing Mg concentration, the layer spacing 
increases significantly for a-V205 (~ 9 % increase from 
XMg = 0 to XMg = 0.5 while using vdW-DF 2 ) but re¬ 
mains fairly constant in 15-V2O5 (~ 2% increase from 
XMg = 0 to XMg = 1 ). However, the layer spacing in 


the 6 phase remains higher than in the a phase across all 
Mg concentrations. Also, the layer spacing seen in the a 
phase (with vdW-DF 2 ) benchmarks better with exper¬ 
imental layer spacing values at low Mg concentrations 
(up to XMg = 0 . 2 ) compared to the 5 phase. Though 
including the van der Waals corrections in DFT leads 
to better agreement with the experimental V2O5 layer 
spacing, the Mg insertion voltage is overestimated^ (by 
18 % as compared to 6% with PBE-I-I/), showing that 
PBE-|-f 7 describes the energetics more accurately than 
vdW-DF 2 . If the Mg2,V205 hull (Figure [^) were to be 
calculated with vdW-DF 2 , we speculate that the energies 
of the demagnesiated structures will shift to higher val¬ 
ues than PBE-|-f 7 , since van der Waals corrections tend 
to penalize under-binded (demagnesiated) structures. 


C. Mg diffusion barriers in V2O5 


Path Distance (%) 



FIG. 4 : (Color online) (a) Activation barriers for Mg diffu¬ 
sion in select limiting cases in a-V205 and (b) for Mg diffusion 
in 5-V2O5 calculated through the NEB method. 

To gain insight into the migration behavior of Mg in 
a and <5 polymorphs, the calculated activation barriers 
using the NEB method are plotted in Figure The 
migration energy is plotted along the diffusion path with 
the energies of the end points referenced to zero and the 
total path distance normalized to 100 %. The diffusion 
paths in both a and 5 polymorphs correspond to the 
shortest Mg hop along the a-direction as in Figure 
and respectively and perpendicular to the h-c plane 
in Figure[^. The energy difference between the site with 
the highest energy along the path (the activated state) 
and the end points is the migration barrier. A simple 
random walk model for diffusion would predict that an 
increase in the activation barrier of ~ 60 meV would 
cause a drop in diffusivity by one order of magnitude at 
298 K. 
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Specifically, we have performed four sets of calcula¬ 
tions: dilute Mg concentration (xMg = 0.08) in the a 
phase (blue dots on Figure!^), high Mg concentration 
(xMg = 0.44) in the a phase (orange triangles), dilute Mg 
concentration (xMg = 0.08) in the <5 phase (red diamonds 
on Figure!^) and high Mg concentrations (xMg = 0.92) 
in the (5 phase (green squares). Due to the mechanical 
instability of the a phase at high Mg concentrations, we 
performed NEB calculations in the e phase. Because the e 
phase has a specific Mg ordering, migration to an equiva¬ 
lent site requires two symmetrically equivalent hops. The 
path in the orange triangles of Figure therefore only 
shows one half of the total path. 

The data in Figure illustrates that the barriers in 
the (5 phase (~ 600 — 760 meV) are consistently much 
lower than in the a-phase (~ 975 — 1120 meV), with 
the respective migration energies adopting “valley” and 
“plateau” shapes. Upon addition of Mg the migration 
barriers in a and S both increase. The differences in the 
magnitude of the migration barriers and the shape of the 
migration energies between the a and 6 can be explained 
by considering the changes in the coordination environ¬ 
ment of Mg along the diffusion path. For example, in the 
a phase, Mg migrates between adjacent 8-fold coordi¬ 
nated sites through a shared 3-fold coordinated site (ac¬ 
tivated state), a net 8—>-3—>-8 coordination change, while 
in the S phase Mg migrates between adjacent 6-fold co¬ 
ordinated sites through two 3-fold coordinated sites sep¬ 
arated by a metastable 5-fold coordinated “valley”, a net 
6 —>-5—>-3—>-6 coordination change. Hence, the lower 
barriers of the 6 phase compared to the a phase are likely 
due to the smaller coordination changes and the higher 
layer spacing in S than a as seen in Figure The indi¬ 
cation of superior diffusivity of Mg in 15 -V 2 O 5 motivates 
investigating the intercalation properties of Mg in the <5 
phase further. 


D. Cluster expansion on Mg in 5-V2O5 and 
temperature-composition phase diagram 

Consistent with the data in Figure all Mg-Va ar¬ 
rangements have higher energy than the linear combina¬ 
tion of < 5 -V 205 and ( 5 -MgV 205 , supporting phase sepa¬ 
ration on the S lattice as illustrated in Figure [^, where 
the zero on the energy scale is referenced to the DFT 
calculated end members of the S phase. A total of 97 
Mg-Va configurations, across Mg concentrations are used 
to construct the CE, which encompasses 13 clusters with 
a RMSE of ~ 9 meV/f.u. The CEs Weighted Cross Val¬ 
idation (WCV) score of ~ 12.25 meV/f.u. indicates a 
very good match with the current input set and good 
predictive capability. In Eigurej^ the staircase plot dis¬ 
plays the error in predicting the formation energies of 
different Mg-Va configurations by the CE against their 
respective DFT formation energies. A good CE will have 
lower errors for configurations that are closer to the hull, 
i.e. shorter absolute distance from the ground state hull. 



Mg Concentration Error (meV/f.u.) 


FIG. 5 : (Color online) (a) DFT and Cluster expansion pre¬ 
dicted formation energies are plotted on the vertical scale with 
respect to different Mg concentrations on the horizontal scale, 
(b) The staircase plot indicates the errors in energies encoun¬ 
tered for structures using the cluster expansion (horizontal 
scale) with respect to their respective distances from the hull 
(vertical scale). 


and higher errors for configurations that are further away 
from the hull. The current CE displays errors below 
10 meV/f.u. for most structures whose formation en¬ 
ergies are smaller than 120 meV/f.u. Also, it can be seen 
in Eigurej^ that the structures with the highest errors in 
the formation energy prediction normally have formation 
energies greater than 125 meV/f.u. 


1 0 I-1-1-1-1-1-1-1-1-1-1-r 



Cluster size (A) 


FIG. 6: (Color online) ECI of the clusters vs. their respective 
cluster size are plotted. The insets (a) and (c) display the 
triplet terms and inset (b) shows the quadruplet term with the 
solid blue lines indicating in-plane interactions and the dotted 
blue lines indicating out-of-plane interactions. All insets are 
displayed on the a-b plane. 


The ECIs for the clusters in the CE, normalized by 
their multiplicity and plotted against their respective 
cluster sizes, are displayed in Eigure The size of a 
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given cluster is indicated by its longest dimension; for ex¬ 
ample, in a triplet the cluster size is given by its longest 
pair. Negative pair terms indicate ‘attraction’ (i.e. Mg- 
Mg and Va-Va pairs are favored) and positive pair terms 
indicate ‘repulsion’ (i.e. Mg-Va pairs are favored). The 
figures inside the graph show the triplets and the quadru¬ 
plet used in the current CE with the solid lines indicating 
interactions in the a-b plane and dotted lines indicating 
interactions out of plane {b is the direction perpendicu¬ 
lar to the V 2 O 5 layers). The orange circles indicate Mg 
atoms. The data in Figure illustrates that the most 
dominant (highest absolute ECI value) cluster of the CE 
is a triplet where Mg ions are along the a-b plane (as 
shown in Figure The most dominant pair term is 
attractive and is the longest pair of the most dominant 
triplet. The negative sign of the dominant triplet and 
the dominant pair terms implies that there are 2 possible 
configurations containing Mg which are stabilized: i) all 
three sites are occupied by Mg, and ii) only one of the 
three sites is occupied by Mg, consistent with the sign 
convention adopted in the CE {at = 1 for occupied Mg 
site and Ui = —1 for a vacancy). 



FIG. 7 : (Color online) Mg-V205 intercalation phase dia¬ 
gram for the S phase. The black line indicates the phase 
boundary between the single and two phase regions obtained 
from Monte Carlo simulations of the CE. 


The temperature-concentration phase diagram for Mg 
intercalation into ( 5 -V 2 O 5 is displayed in Figure The 
black line traces the phase boundary between the single 
and two phase regions, obtained from Monte Carlo sim¬ 
ulations with the numerical hysteresis removed by free 
energy integration. Consistent with the 6 hull in Fig¬ 
ure the Mg-V 2 05 is a phase separating system at 
room temperature with extremely low solubilities at ei¬ 
ther ends (< 1%). Note that only the solid (5-phase is 
considered in this phase diagram. In reality, the high 
temperature part of the phase diagram would probably 
form a eutectic since pure V 2 O 5 melts at ~ 954 K.l^ 


V. DISCUSSION 

In this work, we have performed a first-principles in¬ 
vestigation of Mg intercalation into orthorhombic V 2 O 5 . 
Specifically, we investigated the a and S polymorphs us¬ 
ing DFT calculations, evaluating their respective ground 
state hulls, subsequent voltage curves, and their Mg mi¬ 
gration barriers. For the S polymorph, we constructed 
the composition-temperature phase diagram using the 
CE and CMC approach. The theoretical data we have 
collected sheds light not only on the existing experiments 
intercalating Mg into V 2 O 5 , but also provides a practical 
strategy to improve performance. 

From a thorough comparison of the experimental data 
available in the literature to the calculations performed 
in this work, we conclude that by synthesizing V 2 O 5 and 
intercalating Mg (i.e. beginning in the charged state), 
the structure remains in the a phase. For example, in 
the experimental voltage curveJi^^^IIllHlZI the character¬ 
istic plateau followed by a drop at XMg ~ 0.5 compares 
well with the computed voltage curve for the a phase 
(Figure which shows a similar voltage drop corre¬ 
sponding to the e ordering while ( 5 -V 2 O 5 would show no 
such drop. In X-ray diffraction (XRD) data in the liter¬ 
ature on magnesiated V 2 O 5 , no additional peaks which 
would in dicate t he formation of the 5 phase have been 
observed .I22I2SIMI Also, the observed increase in the layer 
spacin^^ is consistent with the computed predictions of 
layer expansion in the a phase until XMg = 0.5 (Figure]^ 
rather than the b phase which has a minimal increase in 
layer spacing from XMg = 0 to XMg = 1. The migration 
barriers for Mg in the a phase are high (^ 975 meV as 
seen in Figure]^), and indeed, reversible Mg insertion 
can be reliably achieved only when the diffusion length 
is greatly reduced (i.e. in thin films and nano-powders) 
and at very low rates (i.e. ^ 0.5 /iA/cm^ by Gershinsky 
et al^^. Magnesiation past the e-phase (xMg ~ 0.5) is 
expected to be difficult as the potential drops thereby 
reducing the driving force for Mg insertion, and the Mg 
migration barrier increases with Mg concentration in a 
(Figure]^). While the driving force to transform from 
a —(5 is small up to XMg 0.5 (as in Figure [^), it 
steeply increases thereafter, leading us to speculate that 
further magnesiation would lead to the formation of a 
fully magnesiated i 5 -MgV 205 on the surface. 

Our thinking on the magnesiation process of V 2 O 5 is 
summarized in Figure The ground state hull in Fig¬ 
ure suggests that under equilibrium conditions the 
Mg insertion mechanism is through a two-phase reaction, 
by nucleation and growth of magnesiated b phase from 
supersaturated o:, rather than through the metastable 
formation of the e phase. These two reaction pathways 
(cycling between 0 and 50% state of charge) are illus¬ 
trated schematically in Figure]^ with the orange squares 
representing Mg atoms. If nucleation and growth of the 
fully magnesiated b phase (i.e. XMg = 1 ) were to occur, 
there would be no inherent upper limit to magnesium 
insertion up to XMg 1. However, the metastable inser- 
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FIG. 8: (Color online) Possible intercalation pathways for 
Mg in V2O5 up to XMg = 0.5. The left half corresponds to the 
equilibrium case where the S phase nucleates and grows in a 
supersaturated a phase, with a well-defined interface between 
the two phases and the right half corresponds to the Mg atoms 
ordering into the metastable e phase and the lack of a well 
defined interface in this case since e and a have the same 
V2O5 layer stacking. 


tion path of Mg in the a phase, which once fully con¬ 
verted to e phase remains at XMg ~ 0.5, is more consis¬ 
tent with experiments. The reason the system follows 
the metastable insertion path through a is that the equi¬ 
librium path (a-V205 to (5-MgV205), requires structural 
rearrangement of the host structure through the trans¬ 
lation of V2O5 layers, which may kinetically be difficult 
once some Mg is inserted and more strongly bonds the 
layers. Also, a nucleation-growth process involves high 
interfacial energies and may lead to low rates. A similar 
metastable solid solution transformation has been pre¬ 
dicted and documented for other thermodynamic phase 
separating systems 

While our calculations, supported by experimental 
data, suggest that the host V2O5 structure remains in 
the a phase upon Mg intercalation, they also suggest 
an approach to substantially improve the electrochem¬ 
ical properties by cycling Mg beginning in the 6 phase. 
Mg in ^-¥205 not only possesses a higher average voltage 
compared to a (~ 120 mV higher as seen in Figure [^), 
but also a significantly better mobility 600 — 760 meV 
compared to ^ 975 — 1120 meV) which accounts for ap¬ 
proximately 5 orders of magnitude improvement in the 
diffusivity at room temperature (Figure]^. Prior com¬ 
putations have reported higher migration barriers in the 
5 phase compared to the a phase in the char ged limit, 
in contrast to our calculations in Figure |4pil which we 
attribute to the authors allowing only Mg and nearby 
oxygen ions to relax in their NEB calculations. In order 
to cycle Mg in the 5 phase, V2O5 must be prepared in the 
fully discharged state ((5-MgiV205), where the 5 phase 
is thermodynamically stable. Fortunately, the synthesis 
of 5-MgV205 is well established in the literature.!^ 

Since at intermediate Mg concentrations the equilib¬ 
rium state is a coexistence between the demagnesiated 
a-phase and the fully magnesiated i5-phase, the S phase 


must remain metastable over a wide Mg concentration 
range to ensure higher capacities. If the (5-phase is 
not metastable, transformation to the a-phase will take 
place. We speculate that the possibility of S phase 
metastability is likely, given that nucleation and growth 
of the a phase requires restructuring of the host lattice, 
and the absence of mechanically unstable Mg configura¬ 
tions (even at XMg = 0) in ^ (Figure]^) in our calcula¬ 
tions. Also, an applied (over)underpotential is required 
to access a metastable (de)insertion path, which can be 
quantified by the difference between the metastable and 
equilibrium voltage curves in Figure [^. For example, to 
avoid the equilibrium path, an applied underpotential of 
~ 800 mV is required to insert Mg and retain the a-V205 
structure, but only ^ 400 mV is required to remove Mg 
and retain the 5-MgV205 structure, which supports the 
possibility of a metastable S phase. 

Assuming the 5-MgV205 phase remains metastable, 
the temperature-composition phase diagram computed 
for Mg in (5-V2O5 using the CE (Figure indicates a 
phase separating behavior with negligible solubility at 
both end members at room temperature. By investigat¬ 
ing the dominant interactions (ECIs) that contribute to 
the CE, we gain some insight into the possible interca¬ 
lation mechanism. The dominant Mg-Va interactions, 
specifically the triplet and the nearest interlayer pair as 
seen in Figure]^ are entirely contained in the a-b plane, 
which indicates that the (5-V2O5 host lattice will contain 
fully magnesiated and fully demagnesiated domains sep¬ 
arated by an interface along an a-b plane. Hence, Mg 
insertion into the 3D 6-V2O5 structure can be effectively 
described by considering the interactions in each 2D a-b 
plane. 
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FIG. 9: (Golor online) Interplay between the dominant pair 
and triplet terms of the CE stabilizing different Mg-Va ar¬ 
rangements. 


Figure illustrates the interplay between these domi¬ 
nant pair and triplet terms which results in the specific 
sequence of Mg-Va configurations in terms of their rela- 
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tive stability. The orange circles indicate Mg atoms, the 
hollow circles the vacancies, and all insets are viewed in 
the a-b plane. Given the sign convention used in the CE 
(cTi = 1 for Mg and at = —1 for Va) and the negative 
sign of the dominant pair and triplet, the formation of 
Mg-Mg and Va-Va pairs are favored while triplets con¬ 
taining one or three Mg atoms are favored. Thus, a fully 
occupied triplet is most stable due to favorable contri¬ 
butions from both the triplet (~ —40 meV) and the two 
longest pair terms —60 meV in total) resulting in a 
net stabilization of ~ —100 meV, while the triplet with 
two Mg atoms forming the shortest pair and a vacancy 
at the apex is least stable due to unfavorable contribu¬ 
tions from both the pairs and the triplet resulting in a 
destabilizing contribution of ~ -1-100 meV. 


The bottom half of Figure illustrates a sample se¬ 
quence in which Mg atoms fill up sites on a given a- 
h plane. The fully magnesiated structure (right inset) 
is highly stabilized due to the presence of fully hlled 
triplets (~ —100 meV/triplet) while the fully demag- 
nesiated structure (right inset) is stabilized to a lesser 
extent (^ —20 meV/triplet). At an intermediate com¬ 
position, the Mg atoms will arrange themselves in such 
a way that the number of fully filled and one-third Hlled 
triplets (^ —40 meV/triplet, depicted in the centre in¬ 
set) is maximized. Since one-third filled triplets stabilize 
a structure more than triplets containing two Mg atoms, 
non-phase separated configurations at low Mg concentra¬ 
tions (xMg < 0.33) will be more stabilized than those at 
high Mg concentrations (xMg > 0.66), as indicated by 
the higher solubilities at lower Mg concentrations in the 
phase diagram shown in Figure at high temperatures. 


Since the occurrence of fully magnesiated and demag- 
nesiated a-b planes is highly stabilized, the intercalation 
of Mg in the 3D ( 5 -V 2 O 5 structure will then progress via 
propagation of fully magnesiated a-b planes along the c- 
axis. With additional applied overpotential, not only can 
the 6 phase be retained, but also a non-equilibrium solid 
solution intercalation pathway in S can be thermodynam¬ 
ically accessible, leading to further improved kinetics.l^ 
An estimate for the additional overpotential required can 
be computed by considering the lowest energy structure 
at XMg = 0.83 in Figure [§r, whose formation energy is 
53 meV/Mg, resulting in an approximate additional over¬ 
potential requirement of ~ 320 mV. Therefore, the net 
overpotential required to access a solid-solution trans¬ 
formation path entirely in the 6 phase upon charge is 
^ 720 mV, which is comparable to the underpotential 
applied (~ 800 mV) to remain in the metastable a phase 
upon discharge. Hence, we suggest that the electrochem¬ 
ical performance of Mg in V 2 O 5 can be improved by be¬ 
ginning cycling in the discharged state, i 5 -MgV 2 05 , with 
the prospect of improved voltage, capacity, and kinetics. 


VI. CONCLUSIONS 

In this work, we have used first-principles calculations 
to perform an in-depth investigation of Mg intercalation 
in the orthorhombic a and S polymorphs of V 2 O 5 to eval¬ 
uate their suitability as high energy density cathode ma¬ 
terials for Mg-ion batteries. Specifically, we computed 
the ground state hulls and the activation energies for Mg 
migration in both polymorphs. For the S polymorph we 
calculated the temperature-composition phase diagram. 
The equilibrium state of Mga;V 205 (0 < XMg < 1) is 
determined to be a two-phase coexistence between the 
fully magnesiated ( 5 -MgV 205 and fully demagnesiated 
a-V 205 phases. NEB calculations indicate that room- 
temperature Mg migration is several orders of magnitude 
faster in the 5 phase (Em ~ 600 — 760 meV) than in the 
a phase (Em ~ 975 — 1120 meV). 

By comparing the calculated voltage curves and 
changes in the layer spacing with intercalation with avail¬ 
able experimental data on Mg insertion in V 2 O 5 , we con¬ 
clude that the a phase likely remains metastable when 
Mg is initially inserted into fully demagnesiated a-V 205 . 
Although the computed a phase migration barriers indi¬ 
cate poor Mg mobility, consistent with reversible Mg in¬ 
tercalation being achievable exclusively at very low rates 
and in small particles, a-V 2 05 is still one of only three 
known cathode materials where reversible cycling of Mg 
is possible at all (along with chevrel MoeSs and layered 
M0O3). 

Therefore, our finding that the ( 5 -V 2 O 5 polymorph dis¬ 
plays vastly superior Mg mobility as well as a modest 
increase in voltage compared to the a phase is especially 
promising, assuming that the 5 -V 2 O 5 host structure can 
remain stable or metastable across a wide Mg concen¬ 
tration range. Fortunately, the 6 polymorph is thermo¬ 
dynamically stable in the fully discharged state and its 
synthesis procedure well known. 

From our first-principles calculations of the forma¬ 
tion energies of several Mg orderings in the ( 5 -V 2 O 5 
host structure and the resulting computed temperature- 
composition phase diagram, we have also gained insight 
into the possible mechanism of Mg intercalation within 
the 5 host structure. At room temperature, Mg displays 
strong phase-separating behavior with negligible solid- 
solution in the end-member phases and favors the forma¬ 
tion of either completely full or empty a-b planes, which 
are perpendicular to the layers formed by the connecting 
VO 5 pyramids, suggesting an intercalation mechanism 
based on nucleation and growth through the propagation 
of an a-b interface along the c-axis. 
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